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Abstract 

Q_i| The HBT radii extracted in pp and pp collisions at SPS and Tevatron 

show a clear correlation with the charged particle rapidity density. We 
propose to explain the correlation using a simple model where the distance 
from the initial hard parton-parton scattering to the hadronization point 
depends on the energy of the partons emitted. Since the particle multi- 
plicity is correlated with the mean energy of the partons produced we can 
explain the experimental observations without invoking scenarios that as- 
sume a thermal fireball. The model has been applied with success to the 
existing experimental data both in the magnitude and the intensity of the 
correlation. As well, the model has been extended to pp collisions at the 
LHC energy of 14 TeV. The possibilities of a better insight into the string 
spatial development using 3D HBT analysis is discussed. 



1 Introduction 

The size of the source created in pp and pp collisions, as measured with momentum 
correlations, increases with the particle multiplicity The correlation 

of the extracted size with the rapidity density of the collisions, from the HBT 
analysis, was sometimes described as evidence for the existence of a "source" with 
a given size. Some alternative explanations have been given invoking long lived 
resonances and multiple parton interactions jSj- 
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Figure 1: a) Distribution of the transverse jet momentum, for different dN c h/dr] 
ranges. The cut at 3 GeV/c applied in the calculations is visible in the figure, 
b) Distribution of the mean transverse jet momentum (py) versus dN c h/dr]. The 
results were obtained with Pythia at y/s = 1.8 TeV. 

In the present work we are investigating whether the observed behavior may 
be understood in terms of more trivial explanations related to the details of the 
hadronization of the partons leading to jets. We know namely, that the point of 
hadronization of a jet and the point of the initial parton-parton hard scattering 
do not coincide. The distance between them is the so called hadronization length 

[Lhadr)- 

Numerical estimates for the time scale of hadronization vary significantly 
|1] > > [13 , but owing to the Lorentz boost to the laboratory frame, they are pro- 
portional to the energy, Lhadr ~ 0(1) E t (|7j). Hence, there is a dependence of 
Lhadr on the energy due to the Lorentz boost. On the other hand the energy spec- 
trum of the emitted jets depends on the charged particle multiplicity of the events 
as shown in Fig^ Hence if we assume that the hadronization occurs at different 
distances from the initial hard scatterings, depending on the energy of the jet, 
we can expect that this effect may simulate an extended hadronic source without 
invoking the presence of a thermalized source of hadrons. We have followed this 
line of thought in the present work. 



2 Simulation 

The simulation comprises three steps: 

1. The simulation of the particle momenta using a standard Pythia event 
generator. Identification of jets and "underlying event". 
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2. Creation of a spatial distribution of particle origins according to our per- 
ception of the hadronization process. 

3. Implementation of the Bose-Einstein effect and creation of the correlation 
functions. 



2.1 Event generation and jets identification 

With Pythia 6.24 |H] we simulate pp collisions with y/s = 1.8 TeV. The PYCELL 
subroutine, that is part of the generator, is used to identify jets. We applied the 
following parameters: 

• pseudo-rapidity range (77): from -2 to 2 

• number of pseudo-rapidity bins: 1200 

• number of bins in azimuthal angle: 1200 

• threshold transverse energy of particles considered: 

• minimum transverse energy of particles that are used as jet seeds: 0.7 GeV 

• minimum jet transverse energy: 3 GeV 



• maximum jet radius R = \J A<p 2 + Arj 2 : 1 

All particles that do not belong to any jet are treated as an "underlying 
event" . Jet axis pj - the direction along which the jet develops - is defined as 

Ptj = J2 Pti W 

i 

= T^pu (2) 
Ptj 

Vj = ^1^1 (3) 
Ptj 

Pj = {Pxj , Pyj , Pzj ) = Ptj (cos , sin (f)j , sinh ^ ) (4) 

The sums run over all particles that make up a jet. p t i, (pi and rji are transverse 
momentum, azimuthal angle and pseudorapidity of the ith particle, respectively. 



2.2 The Source Models 

The events simulated with Pythia are then treated according to our model. 
Namely, the particles identified above as "underlying events" are given a spatial 
origin centered around the initial hard scattering point, while particles within a 
jet are given spatial coordinates of origin according to one of the models described 
below. 
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Tube It is assumed that the hadronization length (lj, the distance from the 
hard scattering) depends linearly on the initial parton energy that we approx- 
imate by p t j (jet total transverse momentum). Thus, lj = fiPtj, where fi is a 
multiplicative factor that represents our lack of theoretical insight into the pro- 
cess of hadronization. For every parton the loci of hadronization along the jet 
axis (xi) is randomized from Gaussian distributions with a mean equal to lj and 
a 0[ = lj/3, preventing negative values. In the transverse direction (with respect 
to the jet axis) the hadronization points are randomized so that the distance to 
the jet axis follows a Gaussian distribution with a variance equal to at and mean 
value of zero. 



Dynamic width The distribution along the jet axis is the same as above while 
the transverse width at depends linearly on the jet transverse energy and more- 
over, it is a function of the position along the jet axis (see Fig|2I) so the distribution 
of hadronization points is: 

_ fcr t max exp ~ (h ~ Xl) ' 2 if a™ x >af° 

&t\ZhPtj) % _min -f -max ^-min ' ^ ' 

\ u t u u t — °t 

I 2 

where cr™ 8 * = ftPtj, cr™ 111 = 0.5 fm and w = m 2 ^max (w is chosen so a t (xi = 0) and 
a t (xi = 21 j) are equal to cr™ m ). 



For both kinds of geometries the hadronization time is equal to x\. 
The positions of the " underlying event" particles are randomized from a single 
Gaussian distribution with variance a^. Their emission time is always equal to 0. 



2.3 Simulation of BE correlations 

Since the generator does not provide for Bose-Einstein correlations they have to be 
introduced. In our simulation we introduce them using the weighting algorithm 
due to Lednicky p. It is applied during the construction of the correlation 
functions. Each particle pair (i, j) is weighted with a probability pij 

C (Q> K ) = N(&K)J2(i,j) Pij , (6) 

Pij = 1 + cos((p i - Pj) ■ (xj - Xj)) . (7) 

where N(Q, K) is the number of pairs in a given bin, Q = Pi — Pj and K = 
(pi+Pj)/2, Xj and are the 4- vectors of the hadronization points and momentum 
in the pair rest frame, respectively. The probability p^j coincides with the formal 
Born probability density of the Bose-Einstein symmetrized 2-particle plane 
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wave. 



Pij = ^(x^X^p^Pj) ^(XijXyjPi.Pj-) , (8) 

*(p»,Pi,Pi,Pj) = — (e^ x ' + ^ x i + e^ Xi+ip «). (9) 

V2 

2.4 Correlation Functions 

The correlation functions are calculated for particles with \rj\ < 1 and p t > 0.1 
GeV, while no such a constraint is imposed in the jet finding procedure. We 
extract the correlation functions with two types of parameterizations. 

• C(Q t , Q ) = 1 + Xexp(^(Q t 2 R t 2 + Q 2 r 2 )) 

• C(Qout, Qside, Qiong) = ^ + ^^Pi^piQlut^out + Qlide^side + Qfong^iong)) 

where: 

• Qt- is the component of the three-momentum difference perpendicular to 
the three-momentum sum 

• Qq: is the difference of energies 

• Qiong, Qout and Q S ide- are the components of 3-momentum difference vec- 
tor in the Longitudinally Co- Moving System (LCMS). Qi ong is parallel to 
beam, Q S ide is perpendicular to beam and total pair momentum, and Q out 
is perpendicular to Q long and Q side (FigJSJ). 

• i?'s: corresponding radii 

• r: dispersion (radius) in the time domain 

When speaking of the correlation function C(Qt) we mean the projection of 
C(Q t , Qq) on the Q t axis for \Qq\ < 200 MeV. By a double Gaussian fit we mean 
a ID fit with the following form of a correlator 

C{Q) = l + \ ie -(— ) +\ 2 e~y— > (10) 

3 Results 

The application of the model described above, in agreement with the intuition, 
shows that the correlation function indeed changes its shape with increasing 
multiplicity. Using the experimental data we have attempted to adjust the pa- 
rameters of the model. We have fixed at = 0.5 fm what corresponds to the 
typical hadronic size. We found that cr;, = 0.4 fm reproduces the experimen- 
tal results at the lowest multiplicities (E735 has measured Rt = 0.62 fm at 
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direction 



Figure 2: a) Schema of the the "dynamic width" jet geometry model. The spread 
in transverse direction with respect to the jet axis depends on the position along 
jet and its magnitude depends on the jet energy b) Definition of Q ou t, Qside and 



< dN c h/dr] >= 6.75). In the frame of our model it implies a non-negligible 
contribution of hard processes in total particle production even at low multiplic- 
ities. This observation is in an agreement with other observations at Tevatron 

D3a- 

However, we were not able to reach compatibility with the experimental values 
of R t at high multiplicities. The increase of the ft parameter causes a decrease 
of the intercept parameter, while the shape of the correlation function stays 
approximately unchanged. In fact, the width of the peak, thus Rt as given by a 
Gaussian fit even decreases with increasing ft. 

Using the out-side-long (OSL) parametrization we have observed that for this 
model Rout grows with ft, while R S ide stays approximately unchanged. This 
finding matches the intuitive representation. Therefore the applied geometry 
with a constant a t limits the growth of R t . We deduce that R S ide must also 
increase with the jet energy. This led us to the "dynamic width" jet geometry. 

Using that model, we have found that we are able to reproduce the experi- 
mental results with ft = 1.0 and f t = 0.6, see FiglHK and|U This is not a unique 
pair of parameters that gives a good agreement with E735 result. Within some 
range we can decrease fi and find such a value of f t such that we still reproduce 
the experimental result (FigEb)- 

We believe that the precision on the determination of f values could be im- 
proved if results of a 3D HBT analysis were available. Namely, the dependence 
of R out and R S ide on event multiplicity is required. We have found that R s ide 
increases together with f t , and R ou t together with ft (see TableJTJ). It means 
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Figure 3: Cross-section through the 2D distribution of the hadronization points in 
the plane perpendicular to the beam. The difference in height between contribu- 
tions from the "underlying event" (the peak for values around 0) and jets (shoul- 
ders) is a phase space effect. In fact, the majority of particles originate from jets. 
The dynamic width jet geometry with ft = 1.0, ft = 0.6 and < dN c h/dr] >= 12.4. 
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(a) (b) 

Figure 4: a) The Q t correlation function from our model compared with the 
correlation function extracted from [2] (< dN c h/dr] >= 12.5) and b) double 
Gaussian fit to it. The dynamic width jet geometry with ft = 1.0, ft = 0.6 and 
< dN ch /dr] >= 12.4. 
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Figure 5: Qt correlation functions for the case of the dynamic width geometry 
with ft = 1.0 and f t = 0.6 a) dN ch /dr] = 3.2, b) dN ch /dr] = 7.3, c) dN ch /dr] = 
17.2, d) dN ch /dr] = 23.0 
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Figure 6: Results of our model compared to E735 result. The left hand plot 
shows the best fit to the data obtained with fi = 1 and ft = 0.6 while the right 
hand plot demonstrates the variation of the results with varying the / values. 



that we can estimate the size of the jet fragmentation volume using the three 
dimensional correlation analysis. 

In our model the "underlying event" (UE) is somewhat overestimated since 
we have imposed the cut on the jets of less than 3 GeV and this part has been 
added to the UE. We have examined how our results change if we reduce UE by 
removing randomly 50% of particles not assigned to jets. We have found that 
the obtained radii stay unchanged within 10%. 

From FigJ7|we see that the correlation functions extend up to large Q's (0.4 — 
1.0 GeV), similarly to the ones obtained by E735. The slope of the "tail" depends 
on the multiplicity. The simulated correlation function cannot be well represented 
by a single Gaussian as expected from the distribution of particle hadronization 
points shown in Fig|21 As it can be seen e.g. in FigJSt, a better fit of the 
correlation function is obtained using a double Gaussian, albeit it is not yet 
probably the exact form of a correlator for this kind of source. The two radii may 
be understood in term of the smaller one representing the correlations among 
the particles from the "underlying event" and the larger one representing the 
correlations of jet particles with the ones of the "underlying event" , although the 
interplay of the different factors make such a representation only partially true. 

It is important to mention that the extracted radii are much smaller than the 
extent of the source due to the fact that the particles from the "underlying event" 
are traveling while the jet did not yet hadronize! Similarly, particles hadronizing 
first within a jet also moves together with not yet hadronized partons. The 
same argumentation also explains the weak dependence of R t on // in the "tube 
geometry" . 
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Table 1: Dependence of OSL radii on fi and /; for the < dN c h/dr] >= 23.0. The 
fits were made imposing a double Gaussian (Eq. [T0J) on ID projections, taking 
two other components > 50 MeV. In the rows where i?2 is not specified fits did 
not converge and a single Gaussian is used instead. 
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Figure 7: The same as in FigEJ normalized at Qt ~ 500 MeV and fitted with the 
single Gaussian form of the correlator. 
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Table 2: R t radii dependence on dN c h/drj at y/s = 14 TeV, fi = 1.0 and f t = 0.6. 
The correlation functions were fitted with double Gaussian form of correlator (see 
Eq. HOD- 

Fitting such a correlation function with a single Gaussian - as was done in 
E735 - brings large uncertainties (FigEJ), because the obtained result is very 
sensitive to the normalization chosen (at which point correlation function crosses 
1). To be able to compare with the experimental results PQ j2] we have nevertheless 
used the single Gaussian fit. The results of the double Gaussian fits are shown 
in Fig. EJandEl 

Finally we have calculated, using the parameters extracted for the Tevatron 
data, the expected correlation of radii with charged particle multiplicities for the 
maximum LHC energy of 14 TeV. In Table |2] we present results of our model 
obtained at LHC energies. 

4 Conclusions 

Using a simple approach which introduces a dependence of the distance of the 
mean hadronization points of a parton on its energy we have been able to re- 
produce very satisfactorily both the dependence of the radii, and the trend of 
the correlation strength lambda with the rapidity density in pp collisions. The 
present results indicate that there is a possibility of an alternative interpreta- 
tion of the results to those presented in [2] and [T2| where the obtained radii are 
interpreted as evidence for the observation of deconfined matter in pp collisions, 
On the other hand the model a posteriori justifies the hadronization scenario 
envisaged because the free parameter / has been found close to unity for the 
range of multiplicities analyzed. We believe that the LHC with its wider range 
of multiplicities in pp collisions offers interesting possibilities to test our model. 
We have therefore presented here the expected variation of the radii in function 
of charged particle multiplicities at the LHC. Similarly, the effect of the parton 
hadronization should be taken into account in the analysis of HBT radii in heavy- 
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ion collisions [TB*] . 
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